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ABSTRACT 

Context. Heavy-elements, i.e. those beyond the iron peak, mostly form via two neutron capture processes: the slow (s-) and the rapid 
(r-) process. Metal-poor stars should contain fewer isotopes that form via the s-process, according to currently accepted theory. It 
has been shown in several investigations that theory and observation do not agree well, raising questions on the validity of either the 
methodology or the theory. 

Aims. We analyse the metal-poor star HD 140283, for which we have a high quality spectrum. We test whether a three-dimensional 
(3D) local thermodynamic equilibrium (LTE) stellar atmosphere and spectrum synthesis code permits a more reliable analysis of the 
iron abundance and barium isotope ratio than a one-dimensional (ID) LTE analysis. 

Methods. Using 3D hydrodynamical model atmospheres, we examine 91 iron lines of varying strength and formation depth. This 
provides us with the star’s rotational speed. With this, we model the barium isotope ratio by exploiting the hyperfine structure of the 
singly ionised 4554 A resonance line, and study the impact of the uncertainties in the stellar parameters. 

Results. The star’s rotational speed was found to be 1.65 ± 0.05 km s“'. Barium isotopes under the 3D paradigm show a dominant 
r-process signature as 77 ± 6 ± 17% = 0.38 ± 0.02 ± 0.06) of barium isotopes form via the r-process, where errors represent the 

assigned random and systematic errors, respectively. We find that 3D LTE fits reproduce iron line profiles better than those in ID, but 
do not provide a unique abundance (within the uncertainties). However, we demonstrate that the isotopic ratio is robust against this 
shortcoming. 

Conclusions. Our barium isotope result agrees well with currently accepted theory regarding the formation of the heavy-elements 
during the early Galaxy. The improved fit to the asymmetric iron line profiles suggests that the current state of 3D LTE modelling 
provides excellent simulations of fluid flows. However, the abundances they provide are not yet self-consistent. This may improve 
with non-local thermodynamic equilibrium considerations and higher resolution models. 

Key words. Stars: individual: HD 140283 - Stars: Population II - Galaxy: halo - Stars: atmospheres - Techniques: spectroscopic - 
Hydrodynamics 


1. Introduction 

Traditional equivalent width analyses of stellar spectra have em¬ 
ployed one-dimensional (ID) model atmospheres and ignored 
the details of the line profile shapes. However, this ignores valu¬ 
able information on the state of the stellar atmosphere, and for 
some investigations it is the line profile that harbours the sci¬ 
entific reward. In this paper we test whether one class of three- 
dimensional (3D) models are currently better able to model iron 
line profiles in a metal-poor subgiant and consider the barium 
isotope ratio they imply. 

The nearby subgiant star HD 140283 is one of the bright¬ 
est metal-poor stars (V = 7.21, Casagrande et al. 2010) and as 
a consequence of the high signal-to-noise {S /N) and resolving 
power that is achievable, it has been the subject of extensive 

* Based on observations taken using the Subaru Telescope High Dis¬ 
persion Spectrograph (HDS), operated by the National Astronomical 
Observatory of Japan. 

** Colour figures are only available through access of the online jour¬ 
nal. 


spectroscopic study. Its age and origin (Bond et al. 2013) and its 
composition have been independently determined in a plethora 
of publications (e.g. Magain & Zhao 1993; Bonifacio & Molaro 
1998; Thoren & Edvardsson 2000; Asplund et al. 2003; Shchuk¬ 
ina et al. 2005b; Roederer 2012; Lind et al. 2012). It is one of 
the most well-studied metal-poor stars. 

There have been numerous analyses of HD 140283’s 
neutron-capture element abundances (e.g. Peterson 2011; Roed¬ 
erer 2012; Siqueira Mello et al. 2012). By measuring the ra¬ 
tio of a heavy element whose origin is primarily of the rapid 
(r-) neutron-capture process with another synthesised primarily 
through the slow (s-) process, it is possible to determine the ratio 
of r- to s-process in a star. A typical example of this technique is 
the measurement of europium and barium ratios, where 94% and 
81% of the solar system abundance can be traced through the r- 
and s-process, respectively (Arlandini et al. 1999). Using theo¬ 
retical nuclear networks such as Arlandini et al. (1999) or Burris 
et al. (2000) allows one to put limits on the r- and s-process ratio. 
Using the Arlandini et al. (1999) nuclear network, the s-process- 
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only [Ba/Eu]* ratio is [Ba/Eu] = +1.13, while the r-process-only 
ratio is [Ba/Eu] = -0.69. In general, the europium and barium 
ratios indicate that HD 140283 has a higher fraction of r-process 
elements than s-process elements (Lambert & Allende Prieto 
2002; Siqueira Mello et al. 2012). In Gallagher et al. (2010) 
we set a lower limit for the [Ba/Eu] ratio, [Ba/Eu] > -0.66, 
as limited quality in the spectrum around the europium 4129 A 
line and uncertainties concerning blends around the europium 
4205 A line did not allow us to make a definitive measurement. 
However, this result suggests that both an r- and s-process regime 
are possible, considering the limits set in Arlandini et al. (1999). 

Most results agree well with currently accepted theory con¬ 
cerning the origin of the heavy elements in the early Galaxy. 
It was first proposed by Truran (1981), who inferred that the r- 
process should be dominant over the s-process at early times, 
explaining the observations and analysis of several metal-poor 
stars (Spite & Spite 1978). This theory has subsequently shown 
to agree nicely with observational constraints on the s-process as 
a function of metallicity (Eran 9 ois et al. 2007). 

Several studies have examined the isotopic fractions of bar¬ 
ium in HD 140283 (Magain 1995; Lambert & Allende Prieto 
2002; Collet et al. 2009a; Gallagher et al. 2010, 2012), by mea¬ 
suring the asymmetry of the 4554 A resonance line profile at¬ 
tributed to the odd isotopes of barium. The hyperfine structure 
(hfs) varies with the ratio, /odd, which in turn reflects the nucle- 
osynthetic origin of the isotopes^. Eor a 100% r-process mixture, 
/odd = 0.46 and for a 100% s-process mixture, /odd = 0.11 (Ar¬ 
landini et al. 1999). The observational studies found a wide range 
of r- and s-process ratios; Magain (1995) and Gallagher et al. 
(2010) found a high s-process fraction with almost no r-process 
material present, which conflicts with the Truran postulate. The 
Lambert & Allende Prieto (2002) /odd value is intermediate be¬ 
tween an r- and s-process solar system mix. The earlier studies 
used ID stellar atmospheres, but later 3D atmospheres capable of 
studying the impact of velocity fields on line profile asymmetries 
were adopted (Cayrel et al. 2007). Collet et al. (2009a) show dif¬ 
ferences between their ID and 3D results for barium. Prom their 
ID measurement, HD 140283 shows an r-process signature but 
their 3D analysis drives the isotope fraction down to an s-process 
value, /odd = 0.15 + 0.12. One would infer from such results that 
this particular method of attaining the r- and s-process mixture 
is sensitive to the atmospheric modelling process. 

Other methods of determining a star’s barium isotope ratio 
have also been employed. Mashonkina et al. (1999) use non¬ 
local thermodynamic equilibrium (NLTE) model atmospheres 
and spectral synthesis to measure barium abundances. It is im¬ 
portant that the method uses a NLTE approach, as very accu¬ 
rate abundance measurements are required, which a local ther¬ 
modynamic equilibrium (LTE) technique cannot replicate to a 
high enough degree. They measure the abundance in high exci¬ 
tation barium lines (i.e. those other than the barium 4554 and 
4934 A resonance lines) and then force the two resonance lines 
to have the same abundance by increasing or decreasing the odd 
isotope fraction, while maintaining the same abundance in the 
subordinate lines. The ratio of the odd-to-even isotopes is then 
calculated and the r- and s-process ratio can be inferred. This 
method finds a good agreement with the [Ba/Eu] abundance de¬ 
terminations for the same star. Unfortunately, determining the 
/odd ratio in this fashion requires that barium has a high enough 
abundance that the weaker subordinate lines are strong enough to 


‘ [X/Y] = log,o(f®),-log.o(™), 

" /odd = [N('35Ba) + A(‘3^Ba)]/iV(Ba) 


be distinguished from the continuum and indeed the noise in the 
spectrum. Therefore, using such a method to analyse HD 140283 
is not possible as the resonance and subordinate lines are very 
weak; the 4554 A line has an equivalent width W - 20.1 mA; 
the 4934 A line’s equivalent width is 13.6 mA; the subordinate 
lines have an equivalent width range 4 < VT (mA) <13. 

In Gallagher et al. (2010) (henceforth PAPERl) we exam¬ 
ined the barium isotopic fraction in detail by measuring /odd us¬ 
ing the line profile asymmetry technique. In order to do this, 
it was necessary to analyse a number of iron lines so that the 
star’s macroturbulent broadening could be determined, which 
is a necessary “fudge factor” that is included when working 
with ID LTE spectral synthesis codes, since they do not cal¬ 
culate non-thermal motions. Iron was chosen over other ele¬ 
ments as there are numerous lines to choose and it is free 
from hfs. (In Gallagher et al. (2012) we also examined calcium 
lines for this purpose, which agreed to within Icr with the iron 
line determination.) To calculate the star’s macroturbulence, 92 
weak (10 < VT(mA) < 50) and apparently unblended iron 
lines were selected across the 4110-5650 A wavelength range. 
Each line was modelled using KURUCZ06 model atmospheres 
(http://kurucz.harvard.edu/grids.html), with parame¬ 
ters Teff/logg/[Ee/H] = 5750K/3.7/ -2.5, and synthesised us¬ 
ing the ATLAS (Cottrell & Norris 1978) ID LTE spectral synthe¬ 
sis code. They were fit to the observed spectrum while varying 
the abundance, macroturbulence and wavelength zeropoint, us¬ 
ing code, which is discussed at length in PAPERl and Gar¬ 
cia Perez et al. (2009). Large asymmetries were found in the red 
wing of the iron lines’ residuals, which we speculated were due 
to a fundamental limit of ID spectral modelling; the assumption 
of symmetric line broadening, and the neglect of 3D hydrody¬ 
namic motions of the gas. This was to be expected as it is well es¬ 
tablished that high resolution, high signal-to-noise spectra reveal 
asymmetric line profiles (Gray 1980), even in stars with weak 
lines (like HD 140283), that cannot be replicated by the classi¬ 
cal ID LTE approach to modelling a star’s spectrum. In PAPERl 
we speculated but did not attempt to confirm that the asymme¬ 
tries encountered might be better replicated by employing the 
latest 3D time-dependent model atmospheres and spectral syn¬ 
thesis codes. In addtion, it was questioned whether a 3D treat¬ 
ment of the barium 4554 A line would make any improvement 
to the isotope analysis. However, Collet et al. (2009a) analysed 
the barium isotope fraction in both ID and 3D. They found that 
the higher r-process fraction found under ID, /odd = 0.33 + 0.13, 
was severely reduced under 3D to /odd = 0.15 + 0.12, close to 
the solar s-process-only ratio. This suggests that a 3D treatment 
of the barium line increases the inferred relative abundance of s- 
process-only isotopes and, hence reduces the r-process isotopes. 

In the present work we pursue the issue of line asymme¬ 
try by refitting the same set of iron lines used in PAPERl and 
re-examine the barium isotope ratio, using the 3D LTE spec¬ 
tral synthesis code^ Linfor3D (Steffen et al. 2014), which em¬ 
ploys time-dependent 3D model atmospheres produced with 
CO^BOLD (Ereytag et al. 2012). Elemental abundances are com¬ 
puted as a result of our methods, however, the discrepancies be¬ 
tween classical and the 3D methodologies described in this paper 
have encouraged us to be mindful of these values. 

This work is organised as follows: §2 details the observations 
and reduction of the observed HD 140283 spectrum; §3 and §4 
describe the details of the 3D hydrodynamic modelling process 
and examines various properties of the model atmospheres, re¬ 
spectively; §5 discusses the analysis of the iron lines in detail; §6 

^ Link to the LinforSD manual 
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presents the newly calculated /odd value for barium; and finally 
§7 summarises our results. 

2. Observations 

The data are unchanged from PAPER 1: the stellar and thorium- 
argon calibration spectra were obtained during the commission¬ 
ing of the High Dispersion Spectrograph (HDS) mounted on the 
Subaru Telescope. The stellar spectrum represents the sum of 13 
exposures, taken over two nights (22/07/2001 and 29/07/2001), 
with a total exposure time of 82 minutes. This gives a SIN - 
1100 per 12 mA wide pixel around 4500 A, as measured from 
the scatter in the continuum of the reduced spectrum. The resul¬ 
tant stellar spectrum has a wavelength range 4110-5460 A and 
5520-6860 A. The typical resolution as measured from thorium- 
argon lines is R = A/AA - 95 000. The spectrum was reduced 
using the thorium-argon spectrum to wavelength calibrate the 
stellar spectrum, with typical RMS errors of 1.5 mA (Aoki et al. 
2004). 

3. Modelling HD 140283 

The atmosphere parameters used for modelling HD 140283 were 
set to mirror those in the ID LTE study carried out in PAPER 1, 
as closely as possible. The parameters that were used to con¬ 
struct the ID LTE atmospheres were T^g - 5750K, logg = 3.7, 
[Ee/H] = -2.5. The same logg value was selected for the 3D 
atmospheres. Unlike a ID model atmosphere, in a 3D model 
atmosphere the effective temperature is not a control parame¬ 
ter. Instead, values for entropy are set for the material flowing 
into the computational box through the open lower boundary; a 
higher entropy here will lead to a larger flux, and hence higher 
effective temperature (Caffau 2009). Suitable values of entropy 
were selected so that comparable effective temperatures were 
seen for both the ID and 3D model atmospheres. Opacity tables 
for [Ee/H] = -2.5 were not available at the time of this analysis. 
As a result, we computed pairs of model atmospheres, one for 
[Ee/H] = -2.0 and the other for [Ee/H] = -3.0. 

Table 1 tabulates the twelve parameter sets of the CO^BOLD 
(Ereytag et al. 2012) 3D model atmospheres we investigated: six 
of them have been computed for [Ee/H] = -2.0 and the other 
six for [Ee/H] = -3.0. The first four parameter sets that best de¬ 
scribe HD 140283 were computed for this work. The other eight 
with different temperature and logg values were constructed as 
part of the “Cosmological Impact of the Eirst STars” (CIEIST) 
collaboration (Ludwig et al. 2009), and are used to determine the 
sensitivity to derived parameters. Other than when stated, only 
the two high resolution atmospheres are used. 

Three types of atmosphere were computed for each model 
parameter set: a fully 3D model atmosphere computed with 
CO^BOLD (§3.1), an equivalent average 3D model atmosphere, 
referred to as (3D) (see §3.2), and a ID hydrostatic model atmo¬ 
sphere computed by the Lagrangian ID hydrodynamical (LHD) 
model atmosphere code. The reason for computing LHD models 
is noted in §3.3. 

All these models employ opacities based on the MARCS 
stellar atmosphere package (Gustafsson et al. 2008). In §4 we 
compare the 3D and (3D) atmospheres with ID atmospheres 
generated using the LHD package, while in §5-6 we compare the 
resultant 3D synthesis with ID synthesis computed using the ID 
LTE spectral synthesis code ATLAS^ (Cottrell & Norris 1978), 

The Cottrell & Norris code called ATLAS should not be confused 
with the model atmospheres and code of the same name by Kurucz. 


which we employ with ID LTE Kurucz model atmospheres as 
in PAPERl. It is expected that differences in resultant ID synthe¬ 
ses from both atmospheres would be negligible, as it has been 
shown before that differences between the Kurucz and MARCS 
model atmospheres are immaterial (Bonifacio et al. 2009). In¬ 
deed, results from comparisons of both ID techniques show that 
differences in line strengths produced from either technique for 
a given abundance are small for all iron lines analysed in the 
sample - less than 0.4 mA on average, with a line-to-line scatter 
of 0.74 mA. These differences suggest that synthesis using ei¬ 
ther ID method would produce lines with equal abundances for 
a given spectral feature. 

3.1. The 3D C(y BOLD atmospheres 

The 3D models that best describe HD 140283 (i.e. T^g = 5750 K, 
logg = 3.7, [Ee/H] = -2.0 & -3.0) were computed for high 
and low spatial resolutions, using the “box-in-a-star” set-up for 
a small proportion of the stellar atmosphere. Both cover the same 
physical region of a star (covering optical depths 10“® ^ to 10^), 
but the high resolution models have twice as many voxels in each 
Cartesian axis (i.e. eight times as many voxels overall) and opac¬ 
ities are sampled at more than twice the number of frequency 
intervals. The high resolution models were only computed for 
the best prescription of HD 140283 because of the length of time 
required to compute them. All models consist of a series of 20 
temporal structures, which we refer to as snapshots. 

The geometrical size of the computational box is approxi¬ 
mately scaled according to the model’s temperature, logg, and 
metallicity by scaling the box size according to the resultant 
pressure scale height at the surface relative to the solar model 
geometrical size, 5.6x5.6x2.3 Mm. However, metallicity effects 
have a small impact on this scaling factor and are ignored so that 
models that share the same temperature and logg but differ in 
metallicity can be compared without the need to rebin the com¬ 
putational box (Caffau 2009). Effects related to stellar spheric¬ 
ity are neglected as a consequence of the size of the computa¬ 
tional box, relative to the size of the star (radius of HD 140283 
~ 2 Rq^). The computational box corresponds to ~ 0.01% of the 
surface of one hemisphere of HD 140283. 

3.2. The average (3D) CCTBOLD atmospheres 

Average (3D) model atmospheres are produced by spatially av¬ 
eraging the thermal structure of the 3D computational box over 
surfaces of equal Rosseland optical depth (trqss)- This was done 
twenty times, once for each of the 20 selected snapshots that 
make up every atmosphere tabulated in Table 1. Details of the 
procedure that produces CO^BOLD atmospheres, like the ones 
used in this study, can be found in Ludwig & Kucinskas (2012); 
Kucinskas et al. (2013). 

3.3. The 1D model atmospheres 

Every 3D model atmosphere described in Table 1 has a counter¬ 
part ID model atmosphere, produced with the LHD code. They 
were produced using the same atmospheric parameters, chemi¬ 
cal compositions, opacities and equations of state as those used 
to compute the CO^BOLD 3D and (3D) model atmospheres. 
Convection is treated using the Mihalas (1978) formulation of 
the mixing length theory (MLT). Three ID atmospheres were 


^ Based upon logg = 3.7 and M, - 0.8 Mq 
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Table 1. The atmosphere parameters used in the present analysis. Temperatures are given to the nearest Kelvin. The Average is the average of 
20 retained snapshots, while the RMS errors of the 20 snapshots are listed to illustrate typical scatter in effective temperature between snapshots. 
X, y, z represent the axes of the simulated hox, where z represents the vertical dimension. The number of opacity bins used in each model is given, 
and the computed time between the first and last selected snapshots is listed in the column labelled A time. 


Model name 

Desired 

reff(K) 

logg 

[Fe/H] 

Average 

Teff(K) 

RMS 

(reft) 

Spatial resolution 
[x,y,z] 

Geometrical box 
size [x, y, z] (Mm) 

Opacity 

bins 

A time 
(hours) 

d3t57g37mm20n02 

5750 

3.7 

-2.0 

5773 

9 

280 X 280 X 300 

36.6 X 36.6 X 22.2 

14 

5.3 

d3t57g37mm30n02 

5750 

3.7 

-3.0 

5786 

10 

280 X 280 X 300 

36.6 X 36.6 X 22.2 

14 

10.6 

d3t57g37mm20n01 

5750 

3.7 

-2.0 

5777 

11 

140 X 140 X 150 

36.7 X 36.7 X 22.2 

6 

26.4 

d3t57g37mm30n01 

5750 

3.7 

-3.0 

5791 

11 

140 X 140 X 150 

36.7 X 36.7 X 22.2 

6 

26.4 

d3t55g35mm20n01 

5500 

3.5 

-2.0 

5502 

14 

140 X 140 X 150 

49.0 X 49.0 X 35.9 

6 

13.0 

d3t55g35mm30n01 

5500 

3.5 

-3.0 

5536 

12 

140 X 140 X 150 

49.0 X 49.0 X 35.9 

6 

56.0 

d3t55g40mm20n01 

5500 

4.0 

-2.0 

5473 

7 

140 X 140 X 150 

20.1 X20.1 X 10.6 

6 

9.4 

d3t55g40mm30n01 

5500 

4.0 

-3.0 

5475 

6 

140 X 140 X 150 

20.3 X 20.3 X 10.6 

6 

16.9 

d3t59g35mm20n01 

5900 

3.5 

-2.0 

5862 

14 

140 X 140 X 150 

59.3 X 59.3 X 38.7 

6 

31.1 

d3t59g35mm30n01 

5900 

3.5 

-3.0 

5873 

14 

140 X 140 X 150 

59.3 X 59.3 X 38.7 

6 

31.1 

d3t59g40mm20n01 

5900 

4.0 

-2.0 

5857 

7 

140 X 140 X 150 

25.8 X 25.8 X 12.5 

6 

8.3 

d3t59g40mm30n01 

5900 

4.0 

-3.0 

5851 

7 

140 X 140 X 150 

25.8 X 25.8 X 12.5 

6 

8.5 



lOQ "^ROSS 

Fig. 1. Temperature structures of: the full 3D atmospheres based on 280x280x300 grids points (orange density plot); the average (3D) atmosphere 
(solid black line); and the equivalent ID LHD temperature structure (dashed black line) for the [Fe/H] = -2.0 (left) and [Fe/H] = -3.0 (right) 
high resolution atmospheres, for a single snapshot in time. Darker regions in the density plot represent areas in the 3D atmospheres with densely 
populated temperature regions for a given optical depth. Large temperature deviations exist between the (3D) and ID atmospheres in both the deep 
layers (logTRoss > 1-0) of the star and in line forming regions at logTRoss < -2.0, while in layers where T ~ Tcff, the (3D) and ID temperature 
structures are quite similar. 


produced for every parameter set-up described in Table 1, for 
<xmlt = 0.5, 1.0 and 2.0. We discuss the effects that the choice in 
the MLT parameter has on the ID atmosphere in §4.1. 


4. Atmospheric model analysis 

In this section we detail notable deviations between the 
CO^BOLD and equivalent ID LHD model atmospheres. We also 
discuss the impact that chemical composition has on certain pa¬ 
rameters that line formation is most sensitive to. Explanations 
discussed here help to understand some properties that the 3D 
synthesis exhibit, relative to the ID synthesis, which are dis¬ 
cussed throughout §5-6. 


4. 1. Temperature structure 

The temperature structure {T, t relation) of a 3D simulation dif¬ 
fers even along parallel lines of sight. This makes it difficult to 
compare it with a classical ID model atmosphere. Using the 
(3D) temperature structures (§3.2) makes the comparison be¬ 
tween ID and 3D simpler, as Fig. 1 illustrates. 

One can see notable differences in temperature (~ 800 K) be¬ 
tween the (3D) and ID atmospheres at logrRoss < -4.0. While 
the ID model temperature profile begins to plateau, the (3D) 
models continue to cool outwards. Asplund et al. (1999) and As- 
plund & Garcia Perez (2001) attribute this cooling effect to the 
inefficient radiative heating of the outermost layers of low metal- 
licity model atmospheres relative to cooling by overshooting. 
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Fig. 2. The temperature structure differences between the (3D) and ID 
model atmospheres for three values of qtmlt- The (T,r) profiles for the 
[Fe/H] = -2.0 and -3.0 (3D) atmospheres have been averaged to pro¬ 
duce the (3D) reference used in generating this plot. The outer atmo¬ 
sphere is not sensitive to values of ttMLT, while the deeper layers of the 
atmosphere are extremely sensitive to this parameter. 


Similarly, notable differences in temperature appear deeper 
in the star at log tross > 1- This can be controlled by the mixing 
length parameter (aMLx) selected for the ID model. In Fig. 1, 
omlt = 2.0. Selection of different values for aMur in the ID 
model are still unable to adequately reproduce the temperature 
structure seen in the (3D) model (Ludwig & Kucinskas 2012); 
Fig. 2. However, as lines do not form in these regions, it is of 
little consequence to the work presented here. We have illus¬ 
trated the effects that aMcr has on the ID temperature structure 
in Fig. 2 relative to an interpolated (3D) model atmosphere for 
[Fe/H] = -2.5®. While the line forming region in the outer at¬ 
mosphere is relatively unchanged for varying the deeper 

regions of the atmosphere are heavily dependent on this value. 

The large temperature differences between the (3D) and ID 
model atmospheres seen in the outer atmosphere could lead to 
large deviations in line strengths of absorption lines that form 
there, which will lead to large abundance corrections between 
the 3D and ID models. As the grid points in the outermost layer 
of a fully 3D model atmosphere do not all have the same physical 
characteristics, grid points in that layer correspond to a range of 
(low) optical depths. Inspection of the models shows that 84% of 
the voxels in that layer have optical depths logTRoss < -6.0, in¬ 
dicating that model data presented (e.g. Fig. 2) at logrRoss S 
-6.0 should be little affected by numerical boundary effects. 
Deeper layers where absorption lines and the continuum form 
will clearly be unaffected. 

4.2. Ionisation fractions 

The lower temperatures in the outer regions of 3D model atmo¬ 
spheres, relative to the equivalent ID model atmospheres, will 
have an effect on the ratio of neutral to singly ionised iron ab¬ 
sorbers. It is well understood from the Saha equation that in an 
atmosphere such as those used here (Teff = 5750, logg = 3.7), 
Fell is the dominant ionisation stage for iron. Saha’s equa- 

® All other plots presented in this paper are produced by interpolat¬ 
ing the [Fe/H] = -2.0 and -3.0 (3D) (T,t) structures together unless 
otherwise stated. 



log xross 


Fig. 3. (Top panel): Comparisons of the (3D) (solid lines) and ID 
(dashed lines) atmosphere’s ionisation fractions for the neutral and first 
two ionisation levels of iron. The [Fe/H] = -2.0 and -3.0 (3D) at¬ 
mospheres have been interpolated to produce this plot. (Middle panel): 
Ionisation fraction of Fei for the full 3D (contour) atmosphere for 10 
snapshots of the [Fe/H] = -2.0 atmosphere. Contour values repre¬ 
sent non-normalised frequencies. Also included are the equivalent (3D) 
(dash-dot), average 3D (solid line) and ID LHD (dashed line) ionisa¬ 
tion fraction of Fei. (Bottom panel): Same as middle panel for the 
[Fe/H] = -3.0 atmosphere. All ionisation fractions presented in the 
figure were calculated under the assumption of LTE. 


tion relates the first ionised-neutral number density ratio, 
and the electron pressure, P^, to terms depending essentially 
only on the temperature, T, the identity of the atomic species, 
Z, via its ionisation potential, Iz, and partition functions, e.g. 

= f{T,Z). The temperature is almost uniform in the outer 
layers (logTRoss < -1.5) of the ID model, and as a result ^Pe 
is to first order constant there. At lower optical depths, the gas 
pressure falls as the density falls, and with it the source of elec¬ 
trons and hence electron pressure. As falls, ^ must rise, but 
as is already nearly 100% and hence N\ has little room for 
adjustment, instead Aq must fall. Hence, for logTRoss < -L5, 
Ao falls at smaller logTRoss in the roughly isothermal ID atmo¬ 
sphere (Fig. 3). 

Fig. 3 illustrates that the 3D and ID iron fractions are almost 
identical between -1.0 < logTRoss 2.0 (the lower two panels 
show that there are slight deviations), but as the outer 3D atmo¬ 
sphere cools more efficiently, ^Pg - f(T,T) is no longer con¬ 
stant, and falls. At lower optical depths, we see from the lower 
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two panels of Fig. 3 that the neutral fraction of iron is steadily 
increasing over the range -5 < log trqss < -1 • Compared with 
the ID model, the neutral fraction is several times higher. The 
singly ionised fraction is lower to compensate, but only by a 
very small factor (< 5%). This has a noticeable secondary im¬ 
pact on the abundances derived for Fe i lines in 3D, but not of 
course for Fen. One effect (which we will see in §5.3.2) is that 
the Fe i lines in 3D model atmospheres form much further out 
than in ID atmospheres, but Feii lines are less affected. 

5. Modelling iron lines 

In this section we detail the behaviours shown by the 3D synthe¬ 
sis, relative to the ID synthesis presented in PAPERl. 

5.1. Synthesis with LinforSD 

Synthetic spectral line profiles were computed using the 3D 
spectral synthesis code Linfor3D for 91 of the 92 lines analysed 
in PAPERl. The Fel 5107.45 A line was dropped from the 3D 
analysis owing to a blend with another iron line, redward of the 
line centroid. 

Synthesising a 3D spectrum is very time consuming, there¬ 
fore it is impractical to compute a spectrum for the entire ob¬ 
served wavelength range. For a detailed review of the compu¬ 
tational considerations for synthesis with Linfor3D, see Boni¬ 
facio et al. (2013). As such, the computations were limited to 
a small wavelength range of +0.3 A around each line centroid. 
Computational times varied according to the spatial resolution 
of the CO^BOLD atmospheres, the wavelength sampling (set to 
5 X 10“^ A) and the wavelength range over which the synthesis 
was calculated. 

At the time of this research, the opacities necessary for a 
3D model with metallicity [Fe/H] = -2.5 were not available. 
As such, every 3D synthetic prohle used in the present analysis 
represents an interpolation of spectra computed separately for 
model atmospheres with [Fe/H] = -2.0 & -3.0. The abundance 
is altered by changing the tog gf of each line. For each atmo¬ 
sphere, every line was synthesised for 11 abundances (i.e. 11 
log gf values) over a range of 0.4 dex, in steps of 0.04 dex, for 
line strengths that closely matched the observed line’s equiva¬ 
lent width. As such, the log gf values used to produce lines with 
the [Fe/H] = -3.0 model atmosphere were 1 dex larger than 
those used to construct lines with the [Fe/H] = -2.0 model at¬ 
mosphere so that the resulting synthesis from both atmospheres 
would have approximately the same equivalent width. 

In addition, the line profiles were convolved with a Gaussian 
of FWHM = 3.31 km s ', which is used to represent the instru¬ 
mental broadening of the observed spectrum, as ascertained in 
PAPERl. 

5.2. Externai broadening determination 

Analogous to a ID LTE investigation, it is important that all ex¬ 
ternal broadening parameters are adequately resolved so that the 
isotope ratio determination is not affected. From Lambert & Al- 
lende Prieto (2002) and the work conducted in PAPERl, it is 
clear that barium isotope ratios determined in the manner pre¬ 
scribed are extremely sensitive to external broadening effects, as 
^foddl^vpWHM was found to be -0.51 and -0.7 (km s“')“', re¬ 
spectively. 

We have re-analysed all 91 iron lines to determine the rota¬ 
tional speed, i;sin/, of HD 140283, and assume this to account 


for broadening not already modelled by the 3D atmospheres or 
the instrumental broadening. The iron line sample contains 81 
Fe I and 10 Fe ii lines. Similar to PAPERl, each line was individu¬ 
ally fit using ax^ routine, described in Garcia Perez et al. (2009). 
The code determines the best fit synthetic profile by x^ minimi¬ 
sation^, for varying wavelength shift, iron abundance (A(Fe)) 
and, in this work, vsini. The continuum normalisation is also 
calculated and included when determining the three free param¬ 
eters for each line by measuring small regions of the continuum, 
either side of every iron line. 

We model the rotational profile using v sin i profiles as imple¬ 
mented in Linfor3D (Ludwig 2007), which are convolved to the 
synthetic spectra. In total, a grid of 121 synthetic spectra were 
computed for every iron line in the sample, i.e. for 11 values of 

V sin i and 11 values of A(Fe). u sin i values were computed over 
the range 1.0 < vsini (kms^*) < 6.0 with Arisini = 0.5kms“’. 
The iron abundances in the synthetic grid covered a small abun¬ 
dance range of 0.4 dex in steps AA(Fe) = 0.04 dex. Therefore 
the A(Fe) grid was unique to each line as the differences be¬ 
tween iron abundance across the line sample is larger than the 
range of A(Fe) values that needed to be covered by the synthesis 
for an individual line. 

5.2.1. Result 

It was found that vsini - 1.65 + 0.05 km s ', where the error 
represents the random error associated with the line-to-line scat¬ 
ter. In Fig. 4 we have plotted the relationship between v sin i and 
line strength. It is shown that the scatter in v sin i increases at 
smaller line strength, demonstrating the effect of hnite S /N. It 
also shows that the dispersion is asymmetric about the mean 
L/sini value (represented by the dark line) and the scatter, rel¬ 
ative to the mean (given by the lightly shaded region). 

Further testing of the Fe i and Fe ii lines revealed two con¬ 
trasting values for v sin i: the average v sin i value for the Fe i 
lines was 1.72 + 0.05kms“', while the average r/sini value for 
the Fe ii lines was 1.08 + 0.17 km s“', where errors represent the 
standard error associated with line-to-line scatter in each case. It 
is possible that this behaviour is indicative of NLTE influences. 
Similar behaviours were shown in Ee i lines and lithium under 
LTE in Lind et al. (2013) and in lithium in Steffen et al. (2012). 
When NLTE effects were considered in the line synthesis, Fe i 

V sin i values were reduced towards Fe n v sin i values; the NLTE 
profiles were broader than the LTE profiles. 

The corresponding mean iron abundance is A(Fe) = 4.78 ± 
0.05 ± 0.02, where the errors represent the random error asso¬ 
ciated with line-to-line scatter added in quadrature to the uncer¬ 
tainties associated with temperature and gravity determinations 
of HD 140283 (§6.2.3) and the systematic error associated with 
choice in model atmosphere, but as we discuss in §5.3, there 
are considerable systematic features of the iron abundance which 
need further elaboration. 

5.2.2. Individual atmosphere results 

The synthetic iron line spectra produced from the [Ee/H] = -2.0 
& -3.0 model atmospheres (which are both assembled by aver¬ 
aging over the 20 snapshots that make up each atmosphere) were 
also fit to the observed spectrum. The [Ee/H] = -2.0 model at¬ 
mosphere was found to have an iron abundance of 4.76+0.02 dex 

’ Using the definition= E |^(0, - Mif /crf^, where O,, M, and <t, are 
the observed profile, model profile and error estimation in the observed 
profile, cr = (S /Ny\ respectively. 
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Fig. 4. The behaviour of the measured v sin i from the interpolated spec¬ 
trum, relative to the line strength for the 81 Fe i (black dots) and 10 Fe it 
(red dots) lines. The lightest shaded region shows the scatter relative to 
the mean from the sample (cr); the darker shaded region the standard er¬ 
ror (0.05); and the darkest line represents the mean value (1.65 km s“*) 
used in the analysis. 



and y sin/was found to be 1.89+0.05 km s"*. The [Fe/H] = -3.0 
was found to have an iron abundance of 4.79+0.02 dex and v sin / 
was found to be 1.57+0.05 km s“'. No noteworthy improvements 
were seen between the fits for either atmosphere and no other no¬ 
ticeable effects were detected. 


5.3. 3D-1D abundance differences 

It has been established in previous works that deviations in line 
abundance exist between 3D and ID analyses (Asplund et al. 
1999; Asplund & Garcia Perez 2001; Frebel et al. 2008; Kucin- 
skas et al. 2013). The analysis done here likewise found abun¬ 
dance differences in the 3D analysis, relative to the ID analysis 
in PAPERl. Initially we only examine the 81 Fei lines. In §5.3.3, 
we discuss the 10 Fe ii lines. 


5.3.1. Excitation potential 

Fig. 5 shows that there is a clear Fe i abundance dependence on 
excitation potential. From the top panel it is shown that while 
the Fe i lines synthesised under ID show no trend with excitation 
potential, the Fe i lines synthesised under 3D show do a definite 
trend. In the second panel we show the difference between the 
3D and ID abundances. The trend of the Fei abundances with 
excitation potential is extremely well defined; lines with very 
low excitation potential require large reductions in abundance 
under 3D to maintain the same equivalent width as the observed 
profile. In the bottom panel we have shown that there is no ma¬ 
jor correlation between excitation potential and line strength, so 
the trend in the middle panel is indeed an excitation potential 
effect. This effect has been seen in other published works, such 
as Collet et al. (2009b) who examine carbon, nitrogen, oxygen 
and iron abundances in the metal-poor giant HD 122563. They 
speculate that said effect may be due to the differences in the 3D 
and ID temperature structures, shortcomings in their model at¬ 
mospheres and/ or the departures from LTE in the models they 
present. In the following two sections, we demonstrate that the 



Fig. 5. Top panel. The 3D (filled circles) and ID (crosses) abundances 
for the 91 iron lines analysed under LTE. Red open circles and red 
crosses represent the Eeii lines for the 3D and ID synthesis, respec¬ 
tively. While the ID abundances remain fairly consistent for varying ex¬ 
citation potential, the 3D results show a reduction in abundance at low 
excitation potentials. Middle panel. The apparent excitation potential 
dependence shown by the Ee i lines (filled circles), emphasised by plot¬ 
ting the difference between 3D and ID abundances. The Feii lines (red 
open circles) are not affected by this trend. The clear Fe i outlier seen at 
~ 1.5 eV is caused by a poor fit to the iron line under the ID paradigm 
from PAPERl. Bottom panel. The scatter in line strength with excita¬ 
tion potential for the Ee i (filled circles) and Fe ii (red open circles). This 
also shows that there is no line strength dependence on excitation po¬ 
tential and thus the abundance difference in the middle panel is purely 
an effect of excitation potential. 


effect is caused by differences in the thermal structure between 
the ID and 3D models. 

Figs 1 & 2 illustrate that at log+Ross < “4, temperature 
differences approach 800 K. These differences lead to a dra¬ 
matic decrease in the electron pressure in the 3D atmosphere, 
relative to the equivalent ID atmosphere, as shown in the top 
panel of Fig. 6. At the outermost parts of the atmospheres 
(logTRoss < -4.5), the difference in electron pressure is over 
an order of magnitude. In the middle panel of Fig. 6 we have 
plotted the corresponding relative electron densities of the two 
model atmospheres, calculated from the electron pressures, as¬ 
suming Pe = NekT. The large decrease in the electron number 
density in the 3D model, relative to the ID is related to the higher 
neutral fraction seen at log+Ross < -3 in Fig. 3. They suggest 
that the 3D atmospheres have cooled sufficiently to allow elec¬ 
trons to recombine with the ionised iron, whereas temperatures 
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Fig. 6. The <3D) model electron pressure (top panel), electron density 
(middle panel) and gas pressure (bottom panel), relative to the equiva¬ 
lent ID model. 


are still too high in the equivalent ID models. The bottom panel 
of Fig. 6 illustrates that the relative behaviours of the electron 
pressures are not driven by a rapid decrease in the (3D) gas pres¬ 
sures, relative to the ID, as the (3D) gas pressures remain higher 
than the equivalent ID gas pressure in the outer atmosphere. 

It is of course disconcerting that the abundances derived in 
3D have a significant dependence on excitation potential. Simi¬ 
lar trends were found when model atmospheres with higher and 
lower effective temperatures were selected from Table 1 . Berge- 
mann et al. (2012) show that this trend is removed when the 
analysis is conducted under a (3D) NLTE paradigm. This would 
strongly suggest that this unwelcome finding is a feature of the 
assumption of LTE in the line synthesis, and not merely because 
of a poor choice in effective temperature or a feature of the 3D 
(and also (3D)) models. This also implies that the systematic er¬ 
ror associated with A(Ee) under LTE, presented in §5.2.1, is a 
poor estimate of the true systematic error. The excitation poten¬ 
tial was checked against u sin i and no trend was found, implying 
that V sin i remains unaffected by this effect. To explore this fea¬ 
ture further, we investigated the formation depth of the lines. 

5.3.2. Formation 

In this section we discuss details of neutral iron line formation 
for four specific Eei lines. They were chosen to represent four 
extremes of formation that we see in the 81 Eei line sample: a 
high excitation, weak line; a high excitation, stronger line; a low 
excitation, weak line; and a low excitation, stronger line. The 
lines selected are listed in Table 2. 

The 3D abundances tabulated in column (2) of Table 2 rep¬ 
resent values that best fit the observed profiles, calculated in this 
investigation. Their resultant equivalent widths are tabulated in 
column (7). The ID abundances in column (3) represent values 


that best fit the observed lines found in the;^f^ analysis conducted 
in PAPER 1, such that their equivalent widths are roughly equal 
to Wobs in the table. Their resultant equivalent widths are tab¬ 
ulated in column (8) under the heading VTid.id- The equivalent 
widths tabulated in column (9) under the heading 1 Tid,3d do not 
reproduce the ID abundances listed in column (3) of the table, 
but instead represent the resultant ID equivalent widths from the 
LHD synthesis for the 3D abundances tabulated in column (2). 
We note that relative differences between the synthesis produced 
using LHD and KURUCZ06 atmospheres are immaterial. There¬ 
fore, it is of little importance which ID codes produce the values 
tabulated, as both are virtually identical. Prom columns (2) and 
(3) in Table 2, it is clear that large discrepancies appear between 
the 3D and ID abundances for the two low excitation lines. 

The behaviour of the Pei lines can be further explored by 
examining their contribution functions. In Pig. 7 we have plot¬ 
ted the equivalent width contribution functions over a log Tross 
scale for the lines listed in Table 2. Plotting contribution func¬ 
tions in this way allows for direct comparisons of lines with 
varying line strength as the log TRoss-integral of the contribution 
functions represent the lines’ equivalent widths. Por full and pre¬ 
cise documentation on their calculation, see (Steffen et al. 2014). 

The top left and top right plots show the two high excita¬ 
tion Pei lines, 4132.900A (left) and 4143.410A (right). The 
bottom left and right plots show the low excitation 5166.282 A 
and 4461.653 A lines, respectively. The solid lines represent the 
3D contribution functions, while the dashed lines represent the 
ID contribution functions. The triple dashed-dot lines indicate 
the depth at which half the equivalent width has been achieved, 
which we refer to as the characteristic depth of formation. This is 
a useful value to include as it helps to understand where, on av¬ 
erage, the lines are forming in the star. Both the ID and 3D lines 
in Pig. 7 have been computed for the same abundances so that 
their equivalent widths are equal to those in listed in columns (7) 
and (9) in Table 2. 

Pig. 7 and Table 2 demonstrate that the high excitation lines 
show little discrepancy in line strength and abundance in the ID 
and 3D synthesis, or in characteristic depth of formation, relative 
to the low excitation lines. The principal difference, as foreshad¬ 
owed in §4.2, is that because the 3D models are cooler at lower 
optical depths (compared to the almost isothermal ID models), 
the neutral fraction of iron is several times higher than in the 
ID case (Pig. 3) and so the Pei lines begin to form further out. 
Consequently the 3D contribution functions extend over a larger 
range of logTRoss. towards smaller values, and this in turn shifts 
the characteristic formation depth outwards and drive the line 
strength up. To correct for this, one must decrease the number 
of absorbers for the low excitation lines to achieve the correct 
equivalent width, leading to a much lower indicated iron abun¬ 
dance, relative to when they are constructed under ID. 

While the low excitation lines require very little energy to 
excite the lower levels, which means that even collisions in the 
cool outermost layers may be sufficient, the high excitation lines 
require higher energy collisions, which are found only in the 
deeper layers of the star. The ID atmosphere cools only from 
~ 4800 K to ~ 4500 K over the interval -4.0 < logrRoss ^ 
-1.5. Over this interval, the Boltzmann factor (exp-^), which 
gives the population of for example, a 3 eV level relative to a 
0 eV level, reduces only to 62% of its deeper value. Over this 
depth range, the change in the relative populations of 0 eV and 
~ 3 eV levels is slight in the ID model. In contrast, the (3D) 
model drops to ~ 3800 K at log tross = -4.0, and the relative 
population of a 3 eV level compared to a OeV level drops to just 
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Table 2. Various properties of the four selected lines of study. Columns (2) and (3) tabulate the resultant iron abundance found from the 3D analysis 
conducted in the present work, and the ID analysis from paperi. The excitation potential and logg/ values have also been tabulated. Column (6) 
tabulates the observed equivalent widths of the four lines. Column (7) tabulates the equivalent widths of the 3D synthesis, computed with the 3D 
abundances from column (2) (VTsd.sd)- Column (8) tabulates the ID synthesis computed with the ID abundances tabulated in column (3) (Wid.id)- 
Column (9) tabulates the ID synthesis computed with the abundances from column 2 (ITid.sd)- 


t(A) 

A(Fe)3D 

A(Fe)iD 

T 

logg/ 

IFobs 

1F3d,3D 

IFiD.lD 

1Tid,3D 




(eV) 



(mA) 

(mA) 

(mA) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

4132.900 

4.92 

4.99 

2.845 

-1.005 

16.43 

15.62 

15.57 

14.03 

4143.410 

4.81 

4.83 

3.047 

-0.204 

33.49 

34.44 

34.86 

32.51 

4461.653 

4.25 

5.01 

0.087 

-3.194 

40.54 

42.04 

40.39 

12.66 

5166.282 

4.13 

4.97 

0.000 

-4.123 

10.24 

10.36 

10.11 

1.76 
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Fig. 7. Equivalent width contribution functions for four distinct iron lines for the 3D (solid line) and counterpart ID (dashed line) synthesis. 
Integrals of the curves yield the equivalent width of the Fei line, (top left): A high excitation, weak Fei line, (top right): A high excitation, strong 
Fe I line, (bottom left): A low excitation, weak Fe i line, (bottom right): A low excitation, strong Fe i line. These lines were chosen from our sample 
to illustrate the effect that excitation potential has on formation depths of Fei lines. Both the 3D and ID synthesis were given the same input 
transition probability (log gf) and excitation potential in each case. While the high excitation lines show good agreement in line strength, the low 
excitation lines do not. This helps to explain the behaviour shown in Fig. 5. 


15% of its deeper value. In other words, the ratios of the pop¬ 
ulations of the higher excitation levels, relative to those of low 
excitation levels, climbs steeply as you enter the 3D star. They 
increase by a factor > 6 in the 3D case, but only by a factor of 1.5 
in the roughly isothermal ID case. Consequently, the higher neu¬ 


tral iron fraction seen in the outer layers of the 3D star (Fig. 3) 
results in an increase in the opacity in the low excitation lines in 
the outer layers of the (3D) star, but crucially not in the high ex¬ 
citation lines since the high excitation levels are not significantly 
populated until much deeper regions. Thus the high excitation 
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Table 3. Iron abundances of the resultant interpolated ([Fe/H] = -2.5) 
synthesis from every model atmosphere. The corresponding atmosphere 
details can be found in Table 1 . The first row tabulates the results from 
the high resolution model atmospheres used throughout this section, 
therefore we have separated it from the rest of the atmosphere sample. 
Errors tabulated are the standard error in the scatter of the sample. 


TeS 

log^ 

A(Fe) 

A(Fe i) 

A (Fell) 

5750 

3.7 

4.78 + 0.02 

4.74 + 0.02 

5.10 + 0.03 

5750 

3.7 

4.80 + 0.02 

4.76 + 0.02 

5.12 + 0.03 

5500 

3.5 

4.60 + 0.03 

4.55 + 0.03 

5.01 + 0.03 

5500 

4.0 

4.58 + 0.03 

4.50 + 0.03 

5.24 + 0.03 

5900 

3.5 

4.91+0.02 

4.89 ± 0.02 

4.99 + 0.03 

5900 

4.0 

4.83 ± 0.03 

4.77 + 0.02 

5.25 + 0.03 


3D case resembles the ID case, whereas the low excitation 3D 
lines benefit from the higher neutral fraction in the cooler outer 
regions of the 3D models. 

While we understand that the different behaviour of the low 
and high excitation lines in the ID and 3D models results from 
the temperature profiles, the fact that we derive such a strong de¬ 
pendence of A(Fe) on excitation potential remains unpalatable, 
and surely points to shortcomings in the 3D LTE line synthesis 
presented here. An NLTE calculation for iron lines would almost 
certainly decrease the neutral iron fraction and alter level popu¬ 
lations in the outermost regions, where the low excitation lines 
are forming according to the current 3D prescription (see Berge- 
mann et al. 2012). Without attempting an NLTE analysis in this 
work, we do briefly examine the LTE ionisation balance before 
going on to look at line asymmetries as modelled by the current 
3D LTE code. 

The V sin i values determined for all 91 iron lines were com¬ 
pared with the iron line’s characteristic depth of formation and 
no trends were found. Therefore usin/ remains unaffected by 
this particular 3D effect. 

5.3.3. Fe I vs. Fe ii 

Fig. 5 shows that substantial abundance differences exist be¬ 
tween the 2.5 - 3.5 eV Fei and Feii lines (AA(Fe) ~ 0.3 dex), 
which could indicate an incorrect atmosphere parameter assign¬ 
ment. Through analysis using the other atmospheres listed in 
Table 1, it was found that the ionisation equilibrium could be 
shifted. We have tabulated the abundances found from analysis 
of the entire sample of Fei and Feii lines in Table 3, for each 
atmosphere. 

The table shows that the atmosphere with Tetj = 5900 K 
and logg = 3.5 shows considerable improvement in the ionisa¬ 
tion equilibrium. However, from the ID analysis, atmospheres 
with Teff = 5750K and logg = 3.7 showed the closest Fei 
to Fell ionisation equilibrium. Naively, it might seem that se¬ 
lection of a different atmosphere parameters is warranted un¬ 
der 3D. However, none of the atmospheres analysed here show 
any improvement in the excitation potential dependence demon¬ 
strated in Fig. 5, indicating that iron abundances determined un¬ 
der this 3D LTE paradigm are not currently reliable enough, 
and the ionisation equilibrium cannot be established adequately. 
We therefore maintain that the model chosen to best represent 
HD 140283 is that defined in the ID analysis, with = 5750 K 
and logg = 3.7. 


5.4. Line asymmetry analysis 

One of the main aims of this iron line analysis was to determine 
whether the 3D synthesis shows notable improvements to the 
line profile fitting over the ID synthesis from PAPERl. A large 
shortcoming of the ID synthesis from PAPERl was the large red¬ 
wing asymmetry seen in almost all absorption lines, which we 
speculated may be due to stellar convection. 

Under the ID paradigm, all line broadening associated with 
macroturbulent motions in the star (and indeed the star’s rota¬ 
tion) are ignored during line synthesis and instead are treated 
afterwards by convolving a line profile of determined FWHM 
with the synthesis. In PAPERl we used three distinct line pro¬ 
files: a Gaussian, a radial tangential profile, and a v sin i profile 
that is used to model stellar rotation. All these methods assume 
that the broadening, i.e. the redistribution of line opacity caused 
by these turbulent motions, has a symmetric distribution. The 
asymmetries seen in the residuals of the synthetic fits to the ob¬ 
served spectrum in PAPERl show that this is not the case. 

The 3D paradigm is able to approximate the effect of convec¬ 
tion in absorption lines as it replicates a dynamic atmosphere, 
and it is no longer assumed that the redistribution of line opaci¬ 
ties is symmetric. 

In Fig. 8 we have plotted the four Fei lines discussed in 
§5.3.2. They are presented in the same order as Fig. 7 to help 
with direct comparisons. The weak, low excitation line (bottom 
left panel), is well fit by both paradigms, but there is neverthe¬ 
less a small improvement found with the 3D fit in the red wing. 
The high excitation lines (top left and right panels) only show 
small improvements over the ID profiles. For the strong, low ex¬ 
citation line (bottom right panel), clear improvements to the line 
asymmetry in 3D can be seen redward of the line centre, over 
the equivalent ID line profile, particularly as judged by the line 
residuals. In general, similar improvements were found for the 
entire sample of Fe i lines. 

In Fig. 9 we have plotted averaged, i.e. reduced noise, residu¬ 
als for the 3D (black line) and ID (red line) for the entire iron line 
sample (top panel), Fe i lines (middle panel) and Fe ii lines (bot¬ 
tom panel) for comparison. The plots were produced by averag¬ 
ing the residuals of all the iron line best fits. This helps remove 
the random noise contribution found in each individual residual, 
so that only consistent fitting errors remain. Notable improve¬ 
ments can be seen in the 3D analyses. In particular the asymme¬ 
tries redward of the line centre that dominate the ID analysis are 
severely reduced in the 3D analysis. In general then, it can be 
said that the 3D paradigm helps to improve the overall fit to an 
absorption line, as it is better able to reproduce line asymmetries 
associated with a dynamic stellar atmosphere. This result gives 
us some reason to believe that the 3D models do provide a more 
realistic representation of fluid flows in the stellar atmospheres, 
not withstanding the difficulties discussed in §5.3. 

6. Modelling barium and the isotope ratio 

In this section we detail the analysis of the singly ionised barium 
4554 A resonance line. We begin by describing the processes in¬ 
volved in creating the barium line list. 

6.1. Synthesising the barium line 

Barium has seven principal, stable isotopes. The two lightest 
isotopes form via the so-called p-process and have 

negligible contribution to the overall solar barium abundance. 
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Fig. 8. Best fit 3D (solid black lines) and ID (red dashed lines) line profiles compared to the observations (diamonds) for the four Fei lines listed 
in Table 2. A residual plot is shown below each line profile. The same ordering has been used for the plots in this figure as in Fig. 7, i.e. top panels 
are the high excitation lines and the bottom panels are the high excitation lines. 


The other five isotopes are formed through a series of neutron- 
captures via either the s- orr-process mechanisms. The s-process 
can synthesise all of these isotopes, but shielding by 
prevents synthesis of by the r-process. 

The odd isotopes broaden all barium spectral lines because 
of hyperfine structure (hfs), in particular the 4554 A resonance 
line. The odd isotopes contribute to the spectral region closer to 
the wings of the line. The relative strength of the odd isotope 
lines blueward of the line centre is smaller than those redward 
of the line centre. When the r-process fraction is large, the odd 
isotopes contribution to the line profile is increased, and hence 


the overall line profile becomes shallower, broader and the line 
asymmetry is increased also. At the other extreme, when barium 
isotopes form via the s-process, the wing becomes shallower and 
the core stronger. For weak, unsaturated lines (like the 4554 A 
line in HD 140283), the isotope configuration has no effect on 
the overall equivalent width. 


Barium lines are synthesised using hfs information from 
Wendt et al. (1984) and Villemoes et al. (1993). For detailed in¬ 
formation on the line list, we refer the reader to PAPER 1. 
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Fig. 9. Reduced noise residual plot for all 91 iron lines (top panel), 81 
Fe I lines (middle panel) and 10 Fe ii lines (bottom panel) for the ID (red 
line) and 3D (black line) analyses. Notable improvements can be seen 
in both the line core and red wing asymmetries under the 3D analysis, 
which are seen clearly in the ID analysis. Small differences between the 
Fe I and complete iron line sample exist but notable differences between 
the Fe ii and the Fe i lines can be seen. However, this is most likely a 
result of the number of lines in each sample. 


6.2. Barium isotope ratio 

The barium isotope ratio, like usini from the iron line anal¬ 
ysis, was determined using the highest resolution model at¬ 
mospheres in Table 1 labelled as d3t57g37mm20n02 and 
d3t57g37mm30n02. From the analysis of the iron lines (§5.2), 
i/sini = 1.65+0.05kms *. We convolved this f siniprofile with 
the barium profile. The instrumental profile, Vinst = 3.31 kms \ 
which we model as a Gaussian, was also convolved with the bar¬ 
ium profile. Barium profiles were computed with Linfor3D for 
11 barium abundances over a range -1.54 < A(Ba) < -1.14 
with AA(Ba) = 0.04. A total of 11 isotope ratios were synthe¬ 
sised that spanned the entire physical range of isotope config¬ 
urations, i.e. /odd = 0.11 - 0.46, as set by the nuclear network 
presented in Arlandini et al. (1999), with A/odd = 0.0035. Syn¬ 
thetic barium profiles were, like the iron profiles, the result of 
an interpolation of profiles synthesised using the [Fe/H] = -3.0 
and -2.0 atmospheres, in each case with abundances set as if 
[Fe/H] = -2.5. 

The isotope ratio was determined with a routine modi¬ 
fied from that of Garcia Perez et al. (2009), which fits synthetic 
barium profiles with varying barium abundance and isotope con¬ 
figuration. The wavelength shift is also measured by shifting the 
observed profile so that the best fit can be attained, through the 

minimum. It was found that /odd = 0.38, which indicates a 
contribution by the r-process of 77% (see Gallagher et al. 2010, 
Fig. 1 for the /odd to r-process relationship). 

Fig. 10 shows the best fit 3D LTE synthetic barium line. 
Similar to Fig. 8, Fig. 10 illustrates how the 3D and ID spec¬ 


Fig. 10. The best fit 3D (solid black line - /dd = 0.38) and ID (dashed 
red line - /dd = 0.02 from paperi) fits to the observed Ban 4554 A pro¬ 
file (black diamonds). A residual plot is presented in the bottom panel. 
We have also included a lower isotope ratio fit to the observed pro¬ 
file (dashed-dot line) for /dd = 0.25, i.e. 40%, for the same usini and 
A(Ba) values used in the best fit. 


tral calculations differ. Not only is the inferred /odd higher in 3D 
(/odd = 0.02 + 0.06 in PAPERI, as determined from the 4554 A 
and 4934 A line analysis), but the residuals to the fit are also im¬ 
proved, both in the core and the red wing, as noted by inspection 
of Fig. 10 and the two numbers found in the present work 
(X^ = 4.2) and PAPERI (x^ = 6.6). 

6.2.1. Snapshot selection 

In §3.1 we described that each atmosphere presented in this pa¬ 
per represents an average of 20 3D structures in time - snapshots. 
They were selected from the whole run to represent the statis¬ 
tics of the complete ensemble of snapshots closely. Hence, every 
profile analysed in the previous sections is produced from 40 in¬ 
dividual snapshots: 20 from the [Fe/H] = -2.0 atmosphere and 
20 from the [Fe/H] = -3.0 atmosphere. In this section we check 
the influence of the limited statistics and fit individual snapshot 
profiles to the observed spectrum to determine /odd- The sepa¬ 
ration of the snapshots in time is long enough that they can be 
considered as statistically independent samples of the flow field. 

Table 4 tabulates the values of /odd and A(Ba) for the 20 
individual snapshots for each atmosphere. Synthesis times are 
unimportant here, but snapshot number 1 and 20 represent the 
first and last instance in time, respectively. We have used the 
highest resolution atmospheres that best represent HD 140283, 
i.e. Teff = 5750K and logg = 3.7. In this case, usint is fixed to 
the value found from the iron line analysis done with the indi¬ 
vidual atmospheres (i/sint = 1.89 km s“^ for the [Fe/H] = -2.0 
atmosphere and i/sint = 1.57 km s“' for the [Fe/H] = -3.0 at¬ 
mosphere §5.2.2), but is not calculated for individual snapshots 
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Table 4. Values of /odd calculated from the 20 individual snapshots of the [Fe/H] = -2.0 and [Fe/H] = -3.0 atmospheres. The final row of the 
table lists the averages of each column. The errors reported represent the standard error to the mean. 


Snapshot no. 

A(Ba) 

/odd 

A(Ba) 

/odd 

[Fe/H] = 

-2.0 

[Fe/H] = 

-3.0 

1 

-1.38 

0.28 

-1.46 

0.53 

2 

-1.38 

0.37 

-1.45 

0.40 

3 

-1.38 

0.35 

-1.47 

0.39 

4 

-1.38 

0.33 

-1.45 

0.26 

5 

-1.38 

0.27 

-1.46 

0.35 

6 

-1.38 

0.46 

-1.47 

0.42 

7 

-1.38 

0.48 

-1.46 

0.35 

8 

-1.38 

0.57 

-1.46 

0.47 

9 

-1.36 

0.37 

-1.47 

0.53 

10 

-1.39 

0.54 

-1.45 

0.20 

11 

-1.36 

0.46 

-1.47 

0.28 

12 

-1.38 

0.39 

-1.47 

0.49 

13 

-1.38 

0.36 

-1.46 

0.59 

14 

-1.37 

0.41 

-1.45 

0.35 

15 

-1.39 

0.32 

-1.45 

0.28 

16 

-1.37 

0.40 

-1.48 

0.57 

17 

-1.38 

0.40 

-1.48 

0.53 

18 

-1.39 

0.37 

-1.46 

0.33 

19 

-1.39 

0.41 

-1.45 

0.19 

20 

-1.38 

0.46 

-1.46 

0.29 

Average 

-1.38 + 0.002 

0.40 + 0.02 

-1.46 + 0.002 

0.39 + 0.03 


because of the desire to see how different snapshots on a single 
star (fixed v sin i) vary. In addition, each snapshot represents an 
instance in time for the same star, therefore v sin i is no longer 
a free parameter as it will not vary over time. When the values 
for /odd in the table are averaged so that the effective metallic- 
ity is [Fe/H] = -2.5, the resultant value for /odd agrees with the 
value reported in §6.2 at the OAcr level. The abundance, when 
averaged, is identical to that reported in §6.2.2. If vsini were 
calculated using the iron lines for each separate snapshot, /odd 
would most likely remain constant, given that we see this same 
behaviour between /odd and u sin i replicated throughout §6. 

Table 4 indicates an appreciable variation of /odd from snap¬ 
shot to snapshot, reflecting temporal variations of the line profile 
shape, /odd values encompass almost the entire s- to r-process 
configuration possible, and in some snapshots /odd even lies out¬ 
side the range suggested by current nuclear synthesis models. 
Therefore it is essential to average over many the snapshots to 
represent the real star. Table 4 shows that the uncertainty of /odd 
arising from the small box size is reduced to ~ 0.03 by aver¬ 
aging over 20 snapshots, which is ~ 10% of the range in /odd 
(0.11 to 0.46). The error could be reduced further by computing 
additional snapshots. However, the given uncertainty is likely an 
upper limit as the rotational line broadening was not adjusted for 
each snapshot individually. As we shall see in §6.2.3, a case-by¬ 
case adjustment leads to a fairly robust value of /odd- 

6.2.2. Barium abundance differences 

The barium abundance determined from the best fit synthetic 
profile is A(Ba) = -1.43 ± 0.05, where the error assigned rep¬ 
resents the propagated uncertainties in the atmospheric param¬ 
eters (§6.2.3). However, systematic errors encroach on the bar¬ 
ium abundance determined, which are due to model selection. It 


was found that the abundance determined above shifted +0.04 
dex between the [Fe/H] = -2.0 and -3.0 models. Therefore 
A(Ba) = -1.43 + 0.05 + 0.04 as usint has no impact on A(Ba). 

In PAPERl it was found that A(Ba) = -1.28 + 0.08, mean¬ 
ing that an abundance difference of -0.15 dex exists between the 
3D and ID best fit profiles. The Ban 4554 A line is a resonance 
line. Like the low excitation Fei, lines the 4554 A line forms 
further out in the 3D atmosphere and its contribution function 
extends over a larger range of log tross. relative to the ID atmo¬ 
sphere. It was found that the characteristic depth of formation 
for the barium line was log+Ross ~ -2.2 in 3D, while in ID 
logTRoss-1.5. It is therefore reasonable to assume that simi¬ 

lar abundance effects would affect the barium lines, like we have 
seen in the iron lines, and that a 3D NLTE analysis of the bar¬ 
ium abundance may resolve the issues we outlined in §5.3.1 and 
§5.3.2. However, the goal of this work was not to evaluate stellar 
abundances, but to determine the isotope ratio of barium. 

6.2.3. Random uncertainties 

In Table 5 we have tabulated how /odd and A(Ba) are affected 
by changes in temperature and gravity. This was done by syn¬ 
thesising the barium line using each atmosphere listed in Table 1 
and convolving every barium grid with the corresponding v sin i 
value calculated through analysis of the iron lines, which were 
also synthesised using each atmosphere. 

It is shown that /odd is not very sensitive to uncertainties asso¬ 
ciated with temperature and gravity. In PAPERl we state that the 
uncertainty relating to the temperature and log g of HD 140283 
were ctt = +100K and criogg = +0.1. We adopt these uncer¬ 
tainties here as well. The uncertainty associated with the gravity 
of HD 140283 means that logg = 0.003 and the tempera¬ 
ture uncertainty gives CTf^^j - 0.003. In summary, the resul- 
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Table 5. Atmosphere effects on and A(Ba). r/sini has been deter¬ 
mined for each atmosphere via analysis of the resultant iron line sam¬ 
ple synthesis. Results presented here are derived from the interpolated 
model atmospheres, so that the effective [Fe/H] = -2.5. The first row 
presents the results from the high resolution model atmospheres. 


Teff (K) 

log^ 

V sin i 
(kms“^) 

A(Ba) 

/odd 

5750 

3.7 

1.65 + 0.05 

-1.43 

0.38 

5750 

3.7 

2.09 ± 0.05 

-1.41 

0.35 

5500 

3.5 

2.49 ± 0.04 

-1.63 

0.32 

5500 

4.0 

3.19 + 0.03 

-1.49 

0.33 

5900 

3.5 

0.93 ± 0.05 

-1.43 

0.31 

5900 

4.0 

2.51+0.04 

-1.31 

0.33 


tant uncertainty to /odd, derived from the sensitivity to the at¬ 
mosphere parameters (T, log^), is negligible at ±0.004. Mea¬ 
suring how varying vsini directly affects /odd allowed us to 
determine that 6fodd/(>vsmi - -0.31 (kms“')“*. This implies 
that cry;,jj,i,sini = 0.02, as the scatter in i/sinf was found to be 
+0.05 km s ' (§5.2.1). Therefore, the total random error we as¬ 
sign to /odd, derived from vsini, temperature and logg uncer¬ 
tainties, is ±0.02. 

For the barium abundance, we find that cr^(Ba), 7 ' = 0.05 
and cr, 4 (Ba) iogg = 0.02, meaning that we assign A(Ba) the er¬ 
ror ±0.05. We note that because of the effects outlined in §5.3.1 
and §5.3.2, this error is most likely not a true representation of 
the total assigned error on the barium abundance. 

It would appear that /odd is more sensitive to model resolu¬ 
tion, as judged by the differences in the high and low resolution 
model results tabulated in rows (1) and (2), respectively, but this 
is a systematic effect, which we now move on to discuss. 

6.2.4. Systematic uncertainties 

In Table 5 it is shown that v sin i varies quite considerably when 
atmosphere parameters are changed. We find an opposite trend 
between the atmospheric velocity amplitudes in the 3D models 
and the derived v sin i. This suggests that these changes are the 
result of ysinf compensating for the changes in the broaden¬ 
ing by the fluid motions. While this has no immediate impact on 
/odd, as it remains fairly consistent in Table 5, it is useful to quan¬ 
tify the effect that the atmospheric parameters have on v sin i. By 
error propagation, it is found that the systematic error contribu¬ 
tion by ysin/ on /odd is ±0.12. One could conservatively adopt 
±0.12 as an indicative systematic error for an independently de¬ 
rived V sin i. However, this would most likely over estimate the 
total systematic error to /odd, as one would independently derive 
V sin i values for each atmosphere used to calculate /odd, to adjust 
for differences in how the velocity fields are treated under differ¬ 
ent model atmosphere parameters. Thus the systematic error is 
reduced to the scatter in the results for /odd, tabulated in Table 5. 
Therefore, we adopt the systematic error based upon this, which 
is ±0.02. 

Finally, Table 5 shows that there is a systematic uncertainty 
associated with model resolution, although more data would be 
required to make an accurate uncertainty assessment. Neverthe¬ 
less, we estimate the error on /odd to be 0.03, based on what we 
have presented here. The total systematic error on /odd, derived 
from the effects of snapshot selection (0.03 & 0.02 from Table 4), 
model parameters (0.02) and model resolution (0.03), is ±0.06. 


6.3. Result 

We find that /odd = 0.38 ± 0.02 ± 0.06, which means that barium 
shows an r-process signature of 77 ± 6 ± 17% in HD 140283. 
Errors presented represent the total assigned random (§5.2.1 and 
§6.2.3) and systematic errors (§6.2.4), respectively. This not only 
suggests that barium isotopes in HD 140283 have a predomi¬ 
nantly r-process origin, but also strengthens our speculation in 
PAPERl that the classical ID LTE prescription is not adequate 
for the task of determining isotope ratios, using the line asym¬ 
metry method presented in this paper and in PAPERl. We are 
not suggesting that ID LTE modelling techniques are no longer 
viable tools for stellar atmosphere analyses, but rather we are 
highlighting a limitation of the paradigm. 

6.4. Discussion 

In PAPERl it was speculated that a 3D treatment of the convec¬ 
tion would possibly increase the s-process fraction as found by 
Collet et al. (2009a). Clearly the result found here does not agree 
with their determination of /odd and indeed with that speculation 
made in PAPERl . However, it would seem from their v sin i deter¬ 
mination that the atmosphere chosen to model HD 140283 may 
not describe the fluid flows of the atmosphere well, which could 
be the result of several parameters: the numerical resolution of 
the model is not high enough (we see a similar increase in v sin i 
from our lower resolution models); the temperature of the model 
was too low (leading to smaller convective velocities); too few 
opacity bins were used in the model (changing the overall struc¬ 
ture of the atmosphere). One or all of these would result in a 
larger v sin / value to compensate. 

Additionally, we note that the iron lines examined by Collet 
et al. (2009a) are, in general, stronger on average than those used 
in this investigation. We have four iron lines in common with 
Collet et al. Of these four lines, two were weak lines, and one of 
these was an Ee ii line. However, the two stronger iron lines in the 
sample show a higher than average v sin i value, which is a com¬ 
mon feature seen in the stronger lines, as shown in Eig. 4. While 
two lines are not a basis for strong evidence, this could suggest 
that the trend v sin i has with line strength (seen in Eig. 4) is not 
random. If this is the case, then this would suggest that a bias 
in V sin i exists (when synthesising under a 3D LTE paradigm), 
which will depend on the selection criterion used; the stronger 
the lines synthesised under 3D LTE, the larger the v sin i value. 

Einally, Lind et al. (2013), who measure lithium isotopes 
in HD 140283, report vsini - 2.83 ± 0.03 km s' under LTE, 
comparable to that of Collet et al. (2009a). However, the Lind 
et al. published v sin i value they assign under the LTE paradigm 
had been artificially increased to compensate for reduction in 
the instrumental broadening FWHM. The actual v sin i value is 
V sin i - 2.04 ± 0.03 km s“' (Lind 2013). This value agrees at the 
0.8(T level with our i/sini determination and at the 1.8cr level 
with that of Collet et al. (2009a). 

In order to help understand the Ba isotopic ratios derived 
using ID models, we experimented by fitting pure-s and pure-r 
ID synthetic spectra not to noisy observational data but to two 
noise-free 3D synthetic spectra, likewise calculated for pure-s 
and pure-r isotope mixes. The 3D synthetic spectra were thus 
treated heuristically as observational spectra of infinite S /N, and 
small offsets in the wavelength, strength and broadening of the 
synthetic ID profiles were permitted, as was the case when fit¬ 
ting the ID profiles to real observations. 

Interestingly, for both isotopic mixes the 3D profiles were 
better fit by the ID pure-s profiles. In other words, if the 3D 
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synthetic spectra are reasonable approximations to real stellar 
spectra, then this experiment suggests that the subtle differences 
between pure-s and pure-r real spectra are much smaller than the 
differences between ID and 3D spectra, to the extent that fitting 
ID spectra cannot recover the correct isotopic ratio; such fits are 
dominated by the gross ID vs. 3D profile differences rather than 
the subtle s- vs. r-process profile differences. 

This experiment does not demonstrate that current 3D syn¬ 
thetic spectra are a better representation of real spectra, but if 
they are then it helps us understand why ID modelling of a star 
whose heavy-element ratios are r-process-like can lead to ID- 
derived isotopic signatures which appear s-process-like. It also 
underscores the importance of solving remaining issues with 3D 
modelling, such as the dependence of abundance on excitation 
potential discussed extensively in §5.3.1. 


7. Conclusions 

We find that the barium isotopes in HD 140283 show an r- 
process origin as /odd = 0.38 + 0.02 + 0.06 (77 ± 6 ± 17% r- 
process), where the errors represent random and systematic er¬ 
rors, respectively. This result contradicts ID the result presented 
in PAPERl, where the isotope configuration was shown to have 
a distinctly s-process origin, such that /odd = 0.02 + 0.06. It is 
important to note that the ID /odd value reported from PAPERl is 
the result of inverse-variance-weighting the isotope ratios from 
the 4554 A (0.01 ± 0.06) and 4934 A (0.11 + 0.19) lines, which 
were both measured. The differences between the two results are 
most likely due to the way in which a 3D and ID model atmo¬ 
sphere simulate fluid flows in a star. In a ID model, fluid flow 
is approximated by the fudge factor, macroturbulence, which 
must be independently measured so that /odd can be determined. 
It was established in PAPERl that determining macroturbulence 
can lead to problems when measuring /odd as we reported that 
SfoidlSvFWHM - ”0.7 (kms^')^'. For a 3D atmosphere, macro¬ 
turbulence is redundant though v sin i must be derived instead. 
As a 3D atmosphere is dynamic, it is better able to the replicate 
velocity fields of convective cells. Importantly, these fluid flows 
are not symmetric, which leads to asymmetrically broadened line 
profiles. 

The isotopic composition of barium behaves like an asym¬ 
metric broadening mechanism - the larger the r-process signa¬ 
ture, the more asymmetric the line becomes. Therefore /odd is 
extremely sensitive to external broadening mechanisms, in par¬ 
ticular the macroturbulence parameter in a ID analysis and v sin i 
in 3D. Because macroturbulence is measured using a different 
species of line, often iron lines, the macroturbulence will also be 
under or over estimated, impacting the isotope ratio, and hence 
the value for /odd- In addtion, because lines are broadened sym¬ 
metrically with a ID prescription, asymmetric convective effects 
are inevitably mistaken as an isotope signature, leading to incor¬ 
rect estimates of the r-process contribution of /odd- 

It was also demonstrated in §6-4 that 3D barium line profiles, 
for either an s-process-only or r-process-only isotope mixture, 
are better fit by ID s-process-only synthetic profiles. This may 
help to explain why /odd was determined to have an s-process- 
only isotope mixture in PAPERl, if 3D line profiles better repre¬ 
sent real stellar spectra. 

Our result also contradicts the result published by Collet 
et al. (2009a), who also employ a 3D modelling technique to 
calculate /odd in HD 140283. However, we speculate in §6.4 that 
this apparent difference in isotope configuration could be caused 
by a number of differences between the model atmospheres used 


by them and in the present analysis, or by a deeper routed prob¬ 
lem, caused by inherent inadequacies in the 3D LTE line synthe¬ 
sis, though we stress here that this is purely speculative at this 
moment. Nevertheless, it is evident that v sin i, like the macro¬ 
turbulence in a ID investigation, is extremely important when 
measuring /odd- 

We have also shown that in LTE, our 3D model does not 
provide a unique abundance of iron as derived from Eei lines. 
The abundances as derived from different lines have a large de¬ 
pendence on their excitation potential, caused ultimately by the 
3D atmosphere models’ temperature structure. We strongly sus¬ 
pect that an NLTE line formation treatment will resolve this is¬ 
sue. Currently, robust NLTE calculations for iron are at the limit 
of what is computationally feasible. However, Shchukina et al. 
(2005a) and Bergemann et al. (2012) made attempts to tackle the 
3D NLTE problem for iron in HD 140283 by examining line pro¬ 
files produced with 1.5D NLTE and (3D) NLTE line synthesis, 
respectively. Like us, they find a strong trend of the abundances 
from Ee i lines in excitation potential in LTE; in NLTE the trend 
vanishes - in fact is even overcompensated for - making it very 
plausible that NLTE effects are the culprit of the present short¬ 
comings in the abundance and v sin i results we report. 

In comparison to ID, the fits to the observed iron and bar¬ 
ium profiles were, in general, improved under the 3D paradigm. 
Moreover, we demonstrated that changes to the model atmo¬ 
sphere parameters and resolution have little effect on /odd, mak¬ 
ing our result quite robust. However, further work is in order to 
improve the line formation calculations in the 3D models, and 
to clarify the reasons for the conflicting results on the barium 
isotopic composition in HD 140283 reported in the literature. 
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